Entanglement Distribution Revealed by Macroscopic Observations 
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What can we learn about entanglement between individual particles in macroscopic samples by 
observing only the collective properties of the ensembles? Using only a few experimentally feasible 
collective properties, we establish an entanglement measure between two samples of spin-1/2 parti- 
cles (as representatives of two-dimensional quantum systems). This is a tight lower bound for the 
average entanglement between all pairs of spins in general and is equal to the average entanglement 
for a certain class of systems. We compute the entanglement measures for explicit examples and 
show how to generalize the method to more than two samples and multi-partite entanglement. 



Observation of quantum entanglement between in- 
creasingly larger objects is one of the most promising av- 
enues of experimental quantum physics. Eventually, all 
these developments might lead to a full understanding of 
the simultaneous coexistence of a macroscopic classical 
world and an underlying quantum realm. Macroscopic 
samples typically contain N ~ 10 23 particles. Because 
the system's Hilbert space grows exponentially with the 
number of constituent particles, a complete microscopic 
picture of entanglement in large systems seems to be in 
general intractable. The question arises: What can we 
learn about entanglement between constituent particles 
of macroscopic samples, if only limited experimentally 
accessible knowledge about the samples is available? 

There is a strong motivation in addressing this question 
because of recent experimental progress in creating and 
manipulating entangled states of increasing complexity, 
such as spin-squeezed states of two atomic ensembles [l| . 
In such experiments one typically measures only expecta- 
tion values of collective operators of two separated sam- 
ples. It is known that the two samples of spins can be 
characterized as either entangled or separable by mea- 
suring collective spin operators Furthermore, such 
measurements are shown to be sufficient to determine en- 
tanglement measures of Gaussian states [3( and of a pair 
of particles that are extracted from a totally symmetric 
spin state (invariant under exchange of particles) [J]. It 
appears that collective operators cannot be used to fully 
characterize entanglement in composite systems without 
strong requirements on the symmetry of the state. 

Here we present a general and practical method for 
entanglement detection between two samples of spins. It 
solely employs collective spin properties of the samples 
and works irrespective of the number of spin particles 
constituting the samples and with no assumption about 
the symmetry or mixedness of the state. The method 
is based on an entanglement measure which is a tight 
lower bound for average entanglement between all pairs 
of spins belonging to the two samples. This measure is 
equal to the average entanglement for a certain class of 
systems which need not be totally symmetric. We gen- 
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FIG. 1: Two contiguous and non-overlapping spin subsystems 
A and B, each of which contains a large number n of spins. 
What can we learn about entanglement of a pair (a, 0) of 
spins chosen at random where a G A and f3 6 B, if individual 
spins are experimentally not accessible but only the collective 
properties of the samples A and B? What can we learn about 
entanglement between A and B from such measurements? 



eralize the method to obtain the entanglement measure 
between M separated spin samples based on collective 
measurements. The results apply for any entanglement 
monotone that is a convex measure on the set of density 
matrices (e.g., concurrence [H, negativity 0, three-way 
tangle Q)- 

Consider two separated ensembles A and B of spin-i 
particles (figure [J) . Each of them contains a large num- 
ber n of spins. Because of the large dimensions (d = 2") 
of the samples' Hilbert spaces, the structure of entangle- 
ment between the two samples is considerably more com- 
plex than between two single spins. While there are ex- 
perimentally viable methods for detecting entanglement, 
they still require a large number of parameters to be de- 
termined (proportional to d 2 Q). The problem simplifies 
in situations in which each of the ensembles of n spins 
can be treated as one large (total) spin of length n/2. 
This means that, within the ensembles, the individual 
spin-i particles form symmetrized states (Dicke states). 
Though this reduces the dimension of the Hilbert space 
of A (B) to d = n + 1, entanglement determination is 
still demanding for large n both experimentally and the- 
oretically: analytical solutions exist only for pure states 
in general and for mixed states only for small n Q. 

In this paper, we give a method to detect entanglement 
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between large spin samples by measuring only a small 
number of collective spin properties (sample spin compo- 
nents and their correlations), which is independent of the 
sample size n. The collective spin operators are 
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The index i denotes the spatial component of the spins: 
i G {1 = x, 2 = y, 3 = z}. The Pauli matrix of the spin 
at site a £ A is given by a^ a \ and analogously for a 
spin j3 from subsystem B. Note that the collective oper- 
ators satisfy the usual commutation relations [Sf, S^-] = 

\hsijk Sfc , since [af\ay ] ] = 2ie ijh S aa >. 
The spin expectation values and correlations are 
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hij(at,/3), (3) 



and analogously for Sf = \ Ylp^B 9i(fi)> where Sf, Sf G 
\-?f , 2f ], T* B G [-^f-, These are only 15 num- 

bers. Here gi(a), gi{(3) and hij(a,/3) are the (dimension- 
less) expectation values and pair correlations of two single 
spins (a, (3) to which it is assumed there is no experimen- 
tal access: g t {a) = (v\ a) )~, hij(a,/3) = {a[ a) &f ] )~^. 
They are obtained from the actual 4x4 density matrix 
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where 1^ Q ^ is the 2x2 identity matrix in the Hilbert space 
of spin a. 

Out of the experimentally accessible quantities ([2|) and 
(|3|) we will construct a 4x4 density matrix of two virtual 
qubits which describes the collective properties of the two 
spin sets. Its a priori justification is (i) that a general 
treatment of the problem between two large samples of 
spins is intractable because of the high dimensionality, 
and (ii) that we have a fully developed theory of entan- 
glement for two-qubit systems. Therefore, this approach 
is a natural (and successful) way to say something at all 
about the entanglement between two spin systems if only 
collective observables are measured. 

We first introduce the normalized (dimensionless) av- 
erage subsystem expectation values (magnetization per 
particle) and correlations: 
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where sf,tff G [-1,1]. These are the coefficients of the 
virtual density matrix: 
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with a denoting the first and b the second virtual collec- 
tive qubit, associated with subsystems A and B, respec- 
tively. Here, l a , l b , and of are 2x2 identity and Pauli 
matrices for the collective qubits a and b. 

The question is whether the density matrix {?]) is posi- 
tive semi-definite, i.e., whether it is a physical state of two 
qubits. The answer is affirmative and the proof follows 
from consideration of an equal- weight statistical mixture 
of one qubit pair which can be in any of the n 2 states 
Pap. The density matrix of this mixture is the mixture 
of density matrices of all possible pairs (a,/3): p mlx = 
TpTla p Pap- It can easily be seen that /5 m ; x is equal to 
Pab as both are uniquely determined by the same expec- 
tations and correlations: (• 
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Thus, pab is a density matrix. Note that without the nor- 
malizations as given in |5]) and ([6]), the method would not 
work. 

Encapsulated in the following two propositions, we re- 
late the entanglement properties of the virtual qubits to 
those of the spin samples. 

Proposition 1. For any entanglement measure E that 
is convex on the set of density matrices the entangle- 
ment of the virtual density matrix E a b = E(p a b) is a 
lower bound for the average entanglement between all 
pairs E a fj = ^ ^2 a ^ E(p a p). This is an immediate con- 
sequence of the convexity of E: E a b = E(^ J2 a ,/3 Pap) < 

Remarks: First, the result holds for entanglement 
measures that are convex. In certain cases this is di- 
rectly implied by the definition of the entanglement 
measure for mixed states, which involves a convex roof 
E{p) = mm Pii ^ i ^TJi-Pi where the minimiza- 

tion is taken over those probabilities {pi} and pure states 
that realize the density matrix p = X^PilV'iXV'il an d 
E(\ipi) (ipi\) is the entanglement measure of the pure state 
\ipi). Second, the proposition implies that if E a b > 
then at least for one pair we must have E(p a p) > 0. 
Thus, a non-zero value of E a b is a sufficient condition 
for entanglement between the two samples A and B. 
Third, the maximal pairwise concurrence [B[ for symmet- 
ric states is found to be 2/n and is achieved for the W- 
state [ltj. It is conjectured that this remains valid also 
when the symmetry constraint is removed. This suggests 
E a b < E a p < 2/n, if concurrence is used as an entangle- 
ment measure. The existence of this upper bound can be 
seen as a consequence of the monogamy of entanglement. 
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We refer to E a b as pairwise collective entanglement as 
it is determined solely by the expectation and correlation 
values of the collective spin observables. The question 
arises: Under what conditions is E ao equal to the aver- 
age entanglement E a p1 Identifying systems for which the 
equality holds, would allow feasible experimental deter- 
mination of the entanglement distribution in large sam- 
ples by observation of their macroscopic properties only. 
It can easily be seen that, if the state is symmetric under 
exchange of particles within each of the samples, one has 
E a b = E a p = E(p a p) for every pair of particles (a, /?). In 
what follows, we identify an important class of systems 
for which E a b = E a p though the corresponding states 
need not be symmetric under exchange of particles. 

Proposition 2. Consider a system (i) with g z (a) = 
for all a € A and g z (/3) = gf for all [3 € B (translational 
invariance within the subsystems), (ii) with h xx (a,[3) — 
e h yy (a, (3) with e — const = +1 or —1 (all pairs are in 
absolute value equally correlated in the x and y-direction, 
where this correlation may be different in size for differ- 
ent pairs), (hi) with constant sign of the z-correlations, 
i.e., sgn[h zz (a, (3)} = — e for all (a,/?), and (iv) where all 
the remaining expectation values and correlations (g x , 
g y , hij with i ^ j) are zero. Non-vanishing average en- 
tanglement E a p as measured by the negativity [6[ is equal 
to the pairwise collective entanglement E ao , if and only if 
the correlation functions h xx (a, [3) and h yy (a, (3) are each 
constant for all pairs and all pairs have a non-positive 
eigenvalue of their partial transposed density matrix. 

The negativity [6] of a density matrix p is defined as 
E(p) = (Tr|/5 pT | - l)/2, where Tr|/5 pT | stands for the 
trace norm of the partially transposed density matrix 
p pT . Hence the negativity is equal to the modulus of 
the sum of the negative eigenvalues of p pT . 

It is important to stress that proposition 2 holds also 
for states that do not need to be totally symmetric, i.e., 
the h zz (a,f3) may be different for different pairs of par- 
ticles. In general, under the above symmetry, the state 
of the virtual qubit pair is of the form 
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where u±±± = \ (1 ± s^ ± sf ± t a zz) . 
The state of an arbitrary pair (a, (3) of particles has a 
similar structure. For example, the state of an arbitrary 
pair of particles extracted from a spin chain with xxz 
Heisenberg interaction has such a form. The proof of 
proposition 2 is given in the appendix. 

We illustrate the method with some explicit examples. 

1. Dicke states. We consider the Dicke state (general- 
ized W-state) 
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with N > 2 spins, k excitations [1) and N—k non-excited 
spins |0), where < k < N. P$ is the symmetrization 
operator. Within the system we consider two subsystems 
A and B each of size n. Because of the total symmetry 
of the state one has E a j, = E a p = E(p a p) for any size 
ri. Only for the cases where just a single spin (k = 0) or 
all spins are excited (k = N), there is no entanglement 
between two arbitrary pairs or arbitrary sized blocks, re- 
spectively [IH . The global maximum of the entanglement 
E a b, measured by the negativity, is reached for k = ^ 
and its value is E™ x — i jfzj, for all n. It vanishes only 
in the limit N — > oo. 

2. Generalized singlet states. The two subsystems A 
and B, each forming a spin s = are in a generalized 
singlet state: 
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where \m) — |2 s; s + m) denotes the eigenstates of the 
spin operator's z-component. The collective two-qubit 
coefficients are tfP = and the s^ b and t$ with i ^ j 

are all zero. The collective entanglement (negativity) is 
E„h = J- = T- which is non-zero for all sizes n of the 

au 2n 4s 

subsystems and vanishes only in the limit n — » oo. 

3. Generalized singlet state with an admixture of non- 
symmetric correlations. Consider the state 

n 

pMM + (1 -p)® 2 (j 01 U=a<10| + |10> ai/s=a <01 

a— 1 

with p 6 [0,1]. This is a mixture of the generalized singlet 
state (|10|) and n perfectly z-correlated pairs (a, (3 = a). 
This state is not symmetric under particle exchange in 
the zz-correlations. The expectation values s"' b and cor- 
relations tfj with i j remain zero. The correlations 
t xx = t yy are reduced by a factor p compared to those of 
the state (1101). The correlations in z-direction, however, 



are modified and read t"* = — p^V^- — — . Therefore, 
there is a critical number of particles n c = \ ] be- 
yond which there is no collective pairwise entanglement. 



_ 1 l+p—n (1— p) 



> 0. Note 



Only for n < n c we have E ao = 4 
that (jlip is in accordance with proposition 2 (appendix) 
and thus E a b = E a p. Figure [2] shows E a b as a function 
of the spin length s = ^ and the mixing parameter p. 
E a f, is non-zero in regions where p > ff^j and decreases 
inversely proportionally to s. 

Our method can be generalized straightforwardly to 
define multi-partite entanglement of M collective spins 
belonging to M separated samples A\, Am, each con- 
taining a large number of spins n. Any convex multi- 
partite entanglement measure (e.g., M-way tangle) which 
is applied to the corresponding collective matrix of the 
M virtual qubits gives a lower bound for the average 
multi-partite entanglement, obeying the usual constraints 
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FIG. 2: (Color online.) Entanglement E a b between two collec- 
tive spins in the generalized singlet state with an admixture 
of non-symmetric noise as a function of spin length s 
and proportion p of the singlet state in the mixture. The en- 
tanglement is non-zero for sufficiently large p and decreases 
inversely proportionally to s. 



for entanglement sharing such as the Coffman-Kundu- 
Wootters inequality 0]. 

Importantly, in the examples considered above our en- 
tanglement measure scales at most with 1/n and van- 
ishes in the limit of infinitely large subsystem sizes n. 
This is a generic property that follows from the commu- 
tation relation for normalized spins in this limit. Taking 
§i = ±Si = ^J2 a ^i one oh t ams hm„_ >00 [s x , §y] = 
linin^oo i y~ s z = 0. This is sometimes interpreted as 
suggesting that averaged collective observables, like the 
magnetization per particle, represent "macroscopic" or 
classical-like, properties of samples. Note, however, that 
for any n there are n pairs between the subsystems so 
that the number of pairs multiplied by the pairwise col- 
lective entanglement can scale with n, showing the exis- 
tence of entanglement for arbitrarily large n. 

Conclusion. In a recent work it was shown that macro- 
scopic properties such as magnetic susceptibilitycan re- 
veal entanglement within macroscopic samples 12[. The 
present work can be viewed as in a way complementary 
as it demonstrates that macroscopic properties (collec- 
tive spin properties and their correlations) can reveal the 
entanglement distribution between two or more macro- 
scopic samples. On the fundamental side, our method 
demonstrates that there is no reason in principle why 
purely quantum correlations could not have an effect on 
the global properties of objects. On the practical side, it 
enables us to characterize the structure of entanglement 
in large spin systems by performing only a few feasible 
measurements of their collective properties, independent 
of the symmetry and mixedness of the state. 

This work was supported by the Austrian Science 
Foundation, Proj. SFB (No. 1506), the Europ. Commis- 
sion, Proj. QAP (No. 015846), and the British Council 
in Austria. 

Appendix. Proof of proposition 2. Depending on the 
sign e of the zz-correlations, only one eigenvalue of 



and (£T , respectively, can be negative: 
1 
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■ yj(g? + eg*)*+4hl x (a,!3) + e h zz (a, (3) 
V(sS + £^) 2 +4(^) 2 +£t; 
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The corresponding negativities are given by E a p = 
|min(0, fJ, a p)\ and E a b = |min(0, v a b)\- One can express 
v ab as given by 



v a b = M 



A. 



(13) 



where p = Y, a ,p Ma/3 and A = £ Qi/9 [(g*+e£ 

Ahl x (a,p)]^ 2 -l[(s a z +ss b z ) 2 +4(tt x ) 2 ] 1/2 - The quantity 
A is the difference between the entanglement measures 
E a p and E a b, i.e., E a b = E a p — A, for the case that v a b < 
and E a p = 4^ Y, a ,p |min(0, p af3 )\ = |min(0,/2)|. This 
is true, if and only if p a p < for all (a, (3), i.e., all pairs 
are either entangled or have eigenvalue zero. According 
to proposition 1, A is non-negative, i.e., 

sjc 2 + 4 {tf x f < & £ Q)/3 Vc 2 + 4 h 2 xx (a, (3) . (14) 



Here we abbreviated c 
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-es' where the 



latter equal sign is due to |5]| . Inequality GIty becomes an 
equality, i.e., A = 0, if and only if h xx {a,(3) is the same 
for all pairs (a, (3) such that t xx = h xx . Therefore, the 
pairwise collective entanglement E a b equals the average 
entanglement E a p, if and only if for all individual pairs 
fJ-ap < and h xx (a, (3) = e h yy (a, [3) — const for all pairs. 
□ 



[9 
[10 

[11 
[12 



B. Julsgaard, A. Kozhekin, and E. S. Polzik, Nature 
(London) 413, 400 (2001). 

A. S0rensen, L.-M. Duan, J. Cirac, and P. Zoller, Nature 
(London) 409, 63 (2001). 

J. Sherson and K. M0lmer, Phys. Rev. A 71, 033813 
(2005). 

X.Wang and K. M0lmer, Eur. Phys. J. D 18, 385 (2002). 
S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 
(1997). 

K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewen- 
stein, Phys. Rev. A 58, 883 (1998); G. Vidal and R. F. 
Werner, Phys. Rev. A 65, 032314 (2002). 
V. Coffman, J. Kundu, and W. K. Wootters, Phys. Rev. 
A 61, 052306 (2000). 

P. Horodecki and A. Ekert, Phys. Rev. Lett. 89, 127902 
(2002). 

J. Schliemann, Phys. Rev. A 72, 012307 (2005). 

M. Koashi, V. Buzek, and N. Imoto, Phys. Rev. A 62, 

050302(R) (2000). 

V. Vedral, New J. Phys. 6, 102 (2004). 
M. Wiesniak, V. Vedral, and C. Brukner, New J. Phys. 
7, 258 (2005); C. Brukner, V. Vedral, and A. Zeilinger, 
Phys. Rev. A 73, 012110 (2006). 



